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On modeling pressure diffusion 
in non- homogeneous shear flows 

By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 


New models are proposed for the “slow” and “rapid” parts of the pressure diffusive 
transport based on the examination of DNS databases for plane mixing layers and 
wakes. The model for the “slow” part is non-local, but requires the distribution 
of the triple-velocity correlation as a local source. The latter can be computed 
accurately for the normal component from standard gradient diffusion models, but 
such models are inadequate for the cross component. More work is required to 
remedy this situation. 


1. Introduction 

In higher-order turbulence models ‘pressure diffusion’ is usually neglected, or at 
best added to ‘turbulent diffusion’ (Launder 1984) and the two modeled in aggre- 
gate. Pressure diffusion refers to the term 5 , ujp in the Reynolds stress budget; 
turbulent diffusion refers to dkuTujuk- The latter represents the ensemble averaged 
effect of random convection and can often be modeled as a diffusion process; the 
former, however, is harder to explain as diffusion. Turbulent diffusion is usually 
considered to be the dominant diffusion mechanism, and pressure diffusion is con- 
sidered to be negligible. However, Lumley (1975a) showed that for homogeneous 
turbulence the application of symmetry and incompressibility constraints to the ex- 
act equation for the “slow” or non-linear part of the pressure diffusion led to the 
result that its magnitude is 20% that of the triple velocity correlation. In addition, 
it is of opposite sign, so that if turbulent diffusion could be modeled as a gradient 
transport, pressure diffusion would represent counter-gradient transport. Demuren 
et al. (1994) examined DNS databases for several shear flows, namely: the mix- 
ing layer simulation of Rogers and Moser (1994); the wake simulation of Moser et 
al. (1996); the boundary layer simulation of Spalart (1988); and_the backward fac- 
ing step simulation of Le et al. (1993). These confirm for the q 2 -equation that in 
simple shear regions pressure diffusion is roughly 20-30% of turbulent diffusion, and 
it appears to be mostly counter-gradient transport, so that it merely reduces the 
effect of turbulent diffusion, which is mostly gradient transport. Thus, the current 
practice of absorbing pressure diffusion and turbulent diffusion into a single model 
term appears reasonable, as far as the main shear regions are concerned. But the 
DNS data show that near the edges of the shear layers turbulent diffusion decreases 
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FIGURE 1. Reynolds stress budget of q 2 from DNS of a plane wake: , 

production; , turbulent diffusion; , velocity-pressure gradient; , 

pressure diffusion; , viscous dissipation; * , viscous diffusion;* , temporal drift. 

(Inset is a blowup of the vertical axis near the edge). 


rapidly to zero, while pressure diffusion decreases only very gradually, so the latter 
then becomes dominant. Thus, the budgets show that near the free stream edge 
the balance is between pressure transport and mean convection, or temporal drift, 
rather than between turbulent transport and the latter. Further, where shear layer 
interactions occur, as near the middle of a wake, pressure diffusion no longer follows 
counter-gradient transport (see Fig. 1), and Lumley’s model becomes inadequate. 
It in fact becomes additive very close to the center, leading to an overall increase 
in total transport, in contrast to the effect in the simple shear layers on either side. 
Both effects cannot be captured by a mere change in model coefficient. There- 
fore, in order to build a model for total diffusive transport in general shear flows 
one must look beyond the homogeneous model of Lumley. The pressure diffusion 
should also be modeled separately. Inhomogeneous effects can be introduced via 
a non-local model based on the elliptic relaxation concept of Durbin (1991, 1993) 
as implemented in the previous study of Demuren et ai (1994). The local model, 
required in this formulation, will be based on the “slow 11 or non-linear part of the 
pressure diffusion. A splitting of the pressure diffusion into “slow’’ and “rapid 11 
parts is therefore necessary. 

DNS databases for the mixing layer simulation of Rogers and Moser (1994) and 
the wake simulation of Moser et ai (1996) are post- processed to split the velocity- 
pressure gradient and pressure diffusion terms into “slow 11 and “rapid” parts. Sep- 
arate models are then proposed for these terms. 
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JURE 2. Reynolds stress budget of u\ from the DNS of plane wake. (See Fig. 1 
legend.) 
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FIGURE 3. Reynolds stress budget of U 1 U 2 from the DNS of plane wake. (See 
Fig. 1 for legend.) 
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2. Governing equations 

The Reynolds stress equations can be written for time-developing plane shear 
flows as: 


D t UiUj = -( UiU k dkUj + u 3 u k dkUi) - d k u l uju k - (uidj(p/ p) + Ujd t (p/p)) 

— 2vdkUid k Uj + vV 2 UiUj ( 1 ) 

( Dt represents the total derivative and d k the partial derivative in the x k coordi- 
nate.) Thus, the time-derivative is balanced by the production, turbulent transport, 
velocity-pressure gradient correlation, dissipation, and viscous diffusion, respec- 
tively. For the cases under consideration in this study, only normal (£ 2 ) derivatives 
of turbulent statistics are non-zero. There is also only one non- zero mean- velocity 
gradient, fyUi. _ 

Budgets for the u\ normal component of the Reynolds stress and the shear stress 
u\u 2 , obtained from the DNS database for the plane wake of Moser et ai (1996), are 
presented in Figs. 2 and 3, respectively. (In all these figures, y c is the center-line, 8 
is the wake half- width, and 8 m is the mixing layer momentum thickness.) In both 
cases, turbulent and pressure transports are of comparable magnitude. Therefore, 
it would be inconsistent to model one and not the other. And, if one tries to 
model them together, then the comment in the introduction with respect to the in- 
equation also applies to the u^-equation. On the other hand, in the uiu 2 -equation, 
these terms virtually cancel each other out, except near the edges of the shear layers. 
Hence, it appears unlikely that a single composite model could reproduce all these 
features. 

The velocity-pressure gradient correlation can be split into a pressure-strain cor- 
relation and a pressure diffusive transport as: 


-(uidj(p/p) + Ujd,(p/p)) = ( p/p){dju , + d,Uj) - ( l/p){6 j2 d2pu , + 6 j2 d 2 pu } ) (2) 

Further, each of these terms can be split into “slow” and “rapid” parts by splitting 
the pressure p, used in the correlations, into p s and p r , respectively. In the present 
study, p 9 and p r are obtained from solution of the equations: 

V p s — dqUpOpU g “h $22^2 

V 2 p r = — 2d2U\d\u 2 (3) 

Figure 4 shows results of the splittings of the velocity- pressure gradient correla- 
tions into “slow” and “rapid” parts for the four non-zero components of the Reynolds 
stresses. We note that for the diagonal components all energy is produced in the 
streamwise component u\ and then transferred to the normal u\ and transverse 
ul components via pressure scrambling. The transfer mechanism appears to be 
quite different for these components; transfer to the normal component is solely 
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through the “slow” part, and transfer to the transverse component is through the 
“rapid” part. It is conjectured that some structural mechanism must be responsi- 
ble for these, though it could not be identified from the analyzed data. However, 
these results agree quite well with those for the homogeneous shear flow simula- 
tion of Rogers et al (1986). On the other hand, the shear stress results do not 
agree. Whereas the present results for wakes and mixing layers (not shown) show 
the velocity-pressure gradient correlation to be mostly in the “rapid” part, the ho- 
mogeneous shear flow results showed nearly equal distribution between the “slow” 
and “rapid” parts. Figure 5 shows results of the splittings of the pressure diffusive 
transport the u\ and tTfuJ components, all others being zero. In both cases, “slow” 
and “rapid ” parts are comparable, and are mostly of opposite sign. Hence, they 
are, in general, larger in magnitude than the sum. Further, the “slow” part appears 
to represent counter-gradient transport, and the “rapid” part gradient transport, 
i.e., more like turbulent diffusion. 

3. Proposed transport models 

It is proposed to model the transport terms in the Reynolds stress equations 
in three parts, namely, turbulent diffusion, “slow” pressure diffusion, and “rapid” 
pressure diffusion. 

S.l Turbulent diffusion model 

Turbulent diffusive transport (TDIFF) is modeled following the proposal of Mellor 
and Herring (1973), (hereafter denoted MH) a s: 

TDIFF Ui ur — [w iUj u k] == C 3 [( k /£)[(ujMi),j -b(u^U|),j ]j,fc (4) 

(“,*.” represents derivative with respect to £*•) This model is derived from the 
isotropization of coefficients in the more complex Hanjalic and Launder(1972) diffu- 
sion model. It does preserve the symmetry of the indices in the triple- velocity corre- 
lation. It was found by Demuren and Sarkar (1993) to yield the correct anisotropy 
of the Reynolds stress in the wake region of channel flows, and a variant by Mellor 
and Yamada (1986) is widely used in geophysical flows. The model will also be 
used to calculate triple- velocity correlations which are required for the modeling of 
“slow” pressure diffusion in the next section. 

3.2 “Slow” pressure diffusion model 

The “slow” part of the pressure transport (SPDIFF) is modeled using a non-local 
elliptic relaxation approach as: 

SPDIFFujuj = -[tfjfcpTui + SikpJtjU = [6 jk fi + 6 ik fj],k (5) 

where 

L 2 dU, -fi = -f t L (6) 

The local source term f L is given by Lumley’s (1975a) model as: 

ft = P^T, L = ~0.2qhTi 
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(y - yc)/b 

FIGURE 4. Separation of velocity-pressure gradient correlation into “slow” and 

“rapid” parts, from DNS data for plane wake, (a)u^, (b^, (c)?!^, (d)tifu2 : , 

total; , “slow”; , “rapid”. 

It is assumed that the same length scale which governs the non-locality in the 
pressure redistribution would also govern the non-locality in the pressure transport. 
For the u \ u 2 and u\ equations, respectively, q 2 u\ and q 2 ii 2 are obtained from the 
MH model. This treatment represents a generalization of the previous study by 
Demuren et aL (1994) in which the k — e turbulence model was used. 

The principal effect of the non-local model is to “elliptically” spread the influence 
of the local source over the length scale L. Figures 6 and 7 present comparisons of 
the pressure diffusion, based on the local model, to DNS data for mixing layer. In 
each case, two model computations are made; one assumes that the triple-velocity 
correlations are known from DNS, and the other computes them from the MH model. 
For 1 * 2 ? shown in Fig. 6, both computations yield similar results, with peaks that are 
somewhat higher than in the DNS. Hence, full application of the non-local model 
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FIGURE 5. Separation of pressure diffusion into slow and rapid parts: DNS data 
for plane wake, (a )u\, (bJwTwJ. (See Fig. 4 for legend.) 

would produce quite good agreement. For shown in Fig. 7, only the first 

approach, which assumes a pre-knowledge of the triple-velocity correlation, gives 
results in agreement with DNS. The MH model grossly underpredicts the q 2 U\, 
and hence its derivative. This appears to be a general flaw of gradient-diffusion 
models for the triple- velocity correlation. They are usually calibrated to reproduce 
the normal component of the turbulent diffusion in simple shear flows. They fail to 
reproduce other components, if these are present, as in the 3D boundary layer study 
of Schwarz and Bradshaw (1994) or the shearless mixing layer study of Briggs et 
al. (1996). This problem will have to be addressed before a reliable, self-contained 
model for the pressure diffusion of n ] u 2 can be produced. 

3.3 “Rapid” pressure diffusion model 

The “rapid” part of the pressure transport (RPDIFF) is modeled in terms of the 
Reynolds stresses and mean velocity gradients. The simplest such model has the 
form: 


RPDIFF^-uJ = -[AjkPrUi + 6ikPrUj]ik = [Sjk9i + &ik9j]*k (8) 

where 

gi,k = c r uiu~iUi f k ( 9 ) 

Equation (8) is similar in form to the “rapid” part of some pressure-strain mod- 
els. This is consistent with the suggestion of Lumley (1975b) that the traditional 
separation of velocity-pressure gradient correlation into a pressure-strain correla- 
tion and a pressure transport is not unique. Preliminary tests show that c r should 
have a value between 0.1 and 0.3. It has been suggested that this part should also 
be modeled with non-local effects, consistent with the modeling of the slow part 
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'IGURE 6. Model predictions of pressure diffusion of u\ in the plane mixing 

lyer: , DNS; , model with triple-correlations from DNS; , model 

nth triple-correlations from MH. 



( V Vc)l^m 


FIGURE 7. Model predictions of pressure diffusion of u x u 2 in the plane mixing 
layer. (See Fig. 6 for legend.) 
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and the pressure-strain correlation, but evidence for such behavior could not be 
discerned from the DNS data. Further testing is desirable. 

4. Conclusions 

The “slow” and “rapid” parts of the velocity-pressure gradient correlations and 
the pressure transport have been calculated from DNS databases of plane mixing 
layers and wakes. These show that, in agreement with homogeneous shear flow 
simulation, the mechanism for the transfer of energy from the streamwise component 
of the Reynolds stress to the normal component is via the “slow” part, whereas for 
the transverse component it is through the “rapid” part. But pressure transport is 
distributed significantly into both “slow” and “rapid” parts, the former being mostly 
counter- gradient transport, and the latter closer to gradient transport. Models are 
proposed for both parts, which show qualitative agreement with DNS data for the 
normal component of the Reynolds stress but have shortcomings when applied to the 
Reynolds shear stress. The main flaw is the inability of gradient diffusion models 
to predict other than the normal component of the triple-velocity correlation for 
which they have been calibrated. Further development and testing is required. 
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